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1 Introduction 

Our starting point is the Gaussian hypergeometric function F(a, b; c; z) and 
its series representation 



^F(a,b;c;z) = ± ^f\ n z", \z\ < 1, 
r(c) ^ o r(c + n)n\ 

which here is written in terms of Pochhammer symbols 

(x) n = x(x + 1) . . . (x + n — 1) = F(x + n)/T(x). 

The hypergeometric series appears as one solution of the Gaussian (or hyper- 
geometric) differential equation, which is characterized by its three regular 
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singular points at z = 0, 1, oo. The local series solutions at and 1 of this 
differential equation are connected by the continuation formula 



L -F(a, b; c; z) = -J-^ " ^ F(a, b; 1 + a + b - c; 1 - z) 



T(c) v ' ' ' ' T(c-a)T(c-b) 
+ r( r^r^ C) ( 1 - z ) C ' a ~ b nc -a,c-b;l + c-a-b;l-z), (1) 

r(a)r(6) 

(|arg(l-z)| < tt). 

Here we want to show that Eq. ([]]) implies an interesting asymptotic 
expansion for a ratio of products of gamma functions, of which only a special 
case was known before. 

By applying the method of Darboux j|, || to ([[]), we derive in Sec. the 
formula in question. The behaviour of this and a related formula is discussed 
in Sec. and illustrated by a few numerical examples. 



2 Derivation of an asymptotic expansion for 
a ratio of products of gamma functions 

It is well-known that the late coefficients of a Taylor series expansion contain 
information about the nearest singular point of the expanded function 0] . In 
this respect we want to analyze the continuation formula (|IJ), in which then 
only the second, at z — 1 singular term R is relevant, which may be written 

as 

R = T(a + b-c)T(l + c-a-b) ~ (c-a) m (c-fr) m _ „ b+m 
r(o)r(6) ^ T(l + c~a-b + m)m\ [ ' 

By means of the binomial theorem in its hypergeometric-series-form , we may 
expand the power factor 



,,-„-/,+,„ r{a + b-c-m + n j _ „ 

T(a + b — c — m)n 



^c—a—b+m _ 



Z . 



n=0 



Interchanging then the order of the summations and simplifying by means of 
the reflection formula of the gamma function, we arrive at 

1 ^ ^ . (c-a) m (c-b) m T(a + b-c-m + n) n 



This is to be compared with the left-hand side L of ([[]), which is 
T _ 1 f> T(a + n)T(b + n) ^ n 



na)T(b)fo r(c + n)n! 

Comparison of the coefficients of these two power series, which according 
to Darboux |4j] and Schafke and Schmidt [§] should agree asymptotically as 
n — > oo, then yields 

r(q + n)T(b + n) = " (c - q) m (c - 6) _ _ 

T(c + n) i ™ ! 

+ 0(r(a + 6- C- M-1 + n)). 

By means of 

0(T(a + b-c-M-l + n)) = T(a + b-c + n)0{n- M ~ l ) 

and the reflection formula of the gamma function, the relevant formula (0) 
may also be written as 

T(a + n)T(b + n) = i + y (c - a) m (c- b) m ofr^" 1 ) 
T(c + n)T(a + b-c + n) ^ m\(l + c - a - b - n) m 

(3) 

The asymptotic expansion for a ratio of products of gamma functions in this 
form (pj) or the other (0) seems to be new. It is only the special case when 
c = 1 which is known. This special case was stated by Dingle g], first proved 
by Paris0, and reconsidered recently by 01ver|J, who has found a simple 
direct proof. His proof, as well as the proof of Paris, can be adapted easily 
to the more general case when c is different from 1 . Still another proof is 
available || which includes an integral representation of the remainder term. 
Our derivation of Eq. (0) or (Q) is significantly different from all the earlier 
proofs of the case when c = 1. 



3 Discussion and numerical examples 

We now want to discuss our result in the form (|3]). First we observe that the 
substitution c — >• a + b — c leads to the related formula 

T(a + n)T(b + n) = i + V ( a ~ c )^ h ~ c )»* + c^rT^ 1 ) (4) 
r(c + n)T(a + b - c + n) ^ m!(l-c-n) m 
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Which of (||) or (|j) is more advantageous numerically depends on the values 
of the parameters, and in this respect the two formulas complement each 
other. Table [1] shows an example with a set of parameters for which ([$]) gives 
more accurate values than ([|), while Table |2] contains an example for which 
(H) is superior to (§). 

For finite n and M — > oo the series on the right-hand side of (|3|) converges 
if Re(l — c — n) > 0. The same is true for if Re(l + c — a — b — n) > . 
Then, in both cases, the Gaussian summation formula yields 

T(l — c — n)T(l + c — a — b — n) 
n — a — n)T(l — b — n) 

which, by means of the reflection formula of the gamma function, is seen to 
be equal to 

r(a + n)Y{b + n) sin(7r[a + n]) sin(7r[6 + n]) ,^ 

T(c + n)T(a + b — c + n) sin(7r[c + n]) sin(7r[a + b — c + n] ) 

Otherwise @ - @ are divergent asymptotic expansions as n — > oo. 

Although in our derivation n is a sufficiently large positive integer, the 
asymptotic expansions (||) - @j are expected to be valid in a certain sector 
of the complex n -plane, and in fact, the proofs of Paris [^] and of Olver 
apply to complex values of n. 

If the series in (^) or (f|) converge, their sums are equal to (^), which 
generally (if neither c — a nor c — b is equal to an integer ) is different from 
the left-hand side of (§) or (|). Therefore @ and (§) can be valid only in the 
half-planes in which the series do not converge. This means that (^) is an 
asymptotic expansion as n — > oo in the half-plane Re(c— 1+n) > 0, and @) is 
an asymptotic expansion as n — > oo in the half-plane Re(a + 6 — c — 1 + n) > 0. 
Otherwise the series on the right-hand sides represent a different function, 
namely (|5]). 

A few numerical examples may serve for demonstration of these facts. In 
Table [5] , the series converge to (|5|) for n = 10 , and therefore (^) and (f|) are 
not valid. For n = 20, on the other hand, the series diverge and so (^[) and 
(fj) hold. The transition between the two regions is at the line Re(n) = 12.4 
in case of (|3|) or Re(n) = 12.5 in case of In Table [|, we see convergence 
for n = —15 and divergence for n = —5, the transition between the two 
regions being at the line Re(n) = —10.4 in case of (|3]) or Re(ra) = —10.5 in 
case of (g) . 
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M right-hand side of ([$]) right-hand side of(|4]) 



1 


0.9771429 


0.9744681 


2 


0.9773113 


0.9780243 


3 


0.9772978 


0.9769927 


4 


0.9773005 


0.9774980 


5 


0.9772995 «- 


0.9771117 


6 


0.9773001 «- 


0.9775615 


7 


0.9772995 


0.9767519 


8 


0.9773003 


0.9791530 


9 


0.9772983 


0.9652341 


10 


0.9773079 


1.2823765 



exact value of (|) or flf): 0.97729983 

Table 1: Values of the right-hand sides of (§) and (f|) for the parameters 
a = 0.7, 6 = 1.2, c = 0.4. 



M right-hand side of (|j) right-hand side of(f|) 



1 0.968000 


0.972093 


2 0.973760 


0.972350 


3 0.971512 <- 


0.972324 


4 0.973078 


0.972331 


5 0.971231 


0.972327 


6 0.975016 


0.972330 


7 0.959571 


0.972325 


8 1.179434 


0.972342 


9 4.748048 


0.972163 


10 26.430946 


0.968966 


exact value of (0) or (]4 


|): 0.97232850 



Table 2: Values of the right-hand sides of (^) and for the parameters 
a = -0.7, b= -1.2, c= -0.4. 
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M right-hand side of (BI) right-hand side 



n = 10 


1 


0.976000 


0.975000 




2 


0.972434 


0.971912 




3 


0.971341 


0.971037 




4 


0.970882 


0.970687 




5 


0.970651 


0.970517 




6 


0.970520 


0.970423 




7 


0.970440 


0.970367 




8 


0.970388 


0.970331 




9 


0.970352 


0.970307 




10 


0.970326 


0.970290 




exact value of(|3|) 


or (|): 1.94045281 




exact value of 


): 0.97022640 «- 


n = 20 


1 


1.008000 


1.007895 




2 


1.007360 


1.007392 




3 


1.007521 


1.007504 




4 


1.007438 <- 


1.007452 <- 




5 


1.007515 <- 


1.007497 «- 




6 


1.007385 


1.007426 




7 


1.007839 


1.007650 




8 


1.002201 


1.005398 




9 


0.921096 


0.965891 




10 


0.478588 


0.740024 



exact value of (|3|) or 
exact value of (IST): 



d): 1.00747290 
50373645 



Table 3: Values of the right-hand sides of (^) and 
a = -11.7, b = -11.2, c = -11.4 . 



for the parameters 
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M right-hand side of (|3]) right-hand side offQ) 



1 


0.9i 


56667 


0.986957 


2 


0.9i 


55648 


0.985745 


3 


0.9i 


55453 


0.985492 


4 


0.9i 


55397 


0.985415 


5 


0.9i 


55376 


0.985386 


6 


0.9t 


55368 


0.985373 


7 


0.9t 


55363 


0.985367 


8 


0.9^ 


55361 


0.985363 


9 


0.9t 


55360 


0.985361 


10 


0.9i 


55359 


0.985360 



exact value of (|) or (@): 1.97071532 



exact value of (|): 0.98535766 



n = —5 



1 


1.010909 


1.011111 


2 


1.009891 


1.009798 


3 


1.010254 <- 


1.010331 


4 


1.009940 <- 


1.009818 


5 


1.010589 


1.011015 


6 


1.005300 


0.998322 


7 


0.951894 


0.887892 


8 


0.737202 


0.459630 


9 


0.134729 


-0.725230 


10 


-1.243041 


-3.418810 



exact value of (|) or (@): 1.01011438 <- 
exact value of (|): 0.50505719 

Table 4: Values of the right-hand sides of (|3p or for the parameters 
a = 11.7, b = 11.2, c= 11.4. 
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